In-plane hyperbolic polariton tuners in terahertz and long-wave infrared regimes

One of the main bottlenecks in the development of terahertz (THz) and long-wave infrared (LWIR) technologies is the limited intrinsic response of traditional materials. Hyperbolic phonon polaritons (HPhPs) of van der Waals semiconductors couple strongly with THz and LWIR radiation. However, the mismatch of photon − polariton momentum makes far-field excitation of HPhPs challenging. Here, we propose an In-Plane Hyperbolic Polariton Tuner that is based on patterning van der Waals semiconductors, here α-MoO3, into ribbon arrays. We demonstrate that such tuners respond directly to far-field excitation and give rise to LWIR and THz resonances with high quality factors up to 300, which are strongly dependent on in-plane hyperbolic polariton of the patterned α-MoO3. We further show that with this tuner, intensity regulation of reflected and transmitted electromagnetic waves, as well as their wavelength and polarization selection can be achieved. Our results can help the development of THz and LWIR miniaturized devices.

The discovery of two-dimensional (2D) vdW crystals has opened avenues for exploring functional materials and devices in the THz (30-3000 μm) and long-wave infrared (LWIR; 8-15 μm) spectral ranges [1][2][3][4][5][6][7][8][9][10] . The THz and LWIR technologies are of great significance for future photonic and optoelectronic applications, such as 5 G/6 G mobile net-works 11,12 , night vision 13 , biomedical imaging and sensing 14,15 , thermal management 16 , and deep-space exploration 17 . However, their development is always limited by the scarcity of materials with strong and tunable intrinsic optical responses, in particular those used for devices in nanoscale and of room temperature operation.
In the past decades, much effort has been devoted to develop narrow band-gap semiconductors (e.g., mercury cadmium telluride and InSb) and quantum materials (e.g., quantum wells/dots, super lattice) [18][19][20] , whose inter-band, intra-band, or inter-subband optical transitions are found in the LWIR and THz regimes. However, due to the relatively weak light-matter operation, their optical responses are weak, and this is further affected by thermal noise. Thus their devices usually require a cryogenic operation to suppress thermal noise, and the introduction of components such as antennas and/or light absorbing layers to improve the electromagnetic absorption, the dimension of which will be out of the scale of tens of micrometers, not to mention the nanoscale. These will result in complex and largevolume architectures that are not favorable for nanodevices, even miniaturized and portable devices. Moreover, the broadband-and polarization-insensitive optical transitions of these materials also restrict their applications in spectrally and polarization-selective photonic and optoelectronic devices, which are in particular interesting in a modern information society. Metamaterials and metasurfaces comprised of artificially designed metallic or dielectric unit cells are able to confine the THz and LWIR electromagnetic waves to enhance light−matter interactions, which therefore give rise to a variety of functional devices [21][22][23][24][25][26] . However, despite intense research efforts, in these spectral ranges, the unit cells of many conventional metamaterial/metasurface often offer restricted confinement due to high losses 21,27 . These are not conducive to device integration, and will also increase device power consumption.
However, so far, the demonstration of utilizing the above exotic characteristics in a practical device is not given. The challenge lies in that one needs to compensate for the large photon−polariton momentum mismatch for far-field excitation and far-field characterization of the HPhPs. In the previous studies, the HPhPs of vdW crystals have been observed by near-field nano-imaging techniques [8][9][10]28 , relying on using a metallic nanotip to compensate for the large momentum mismatch between free-space photons and polaritons. For most practical device applications, direct excitation of the HPhPs from the far-field is necessary. Some earlier studies indicate that it is possible to pattern the surface of vdW crystals, such as with graphene [40][41][42] , hexagonal boron nitride 7,43,44 , semi-metals 45 , and topological insulators 46 to excite and measure the various types of polaritons. However, these previous studies focused on plasmon polaritons with in-plane isotropic [40][41][42]46 and hyperbolic dispersions 45 , as well as phonon polaritons with in-plane isotropic dispersion 7,43,44 . Far-field excitation and characterization of the tunable in-plane HPhPs in vdW crystals, especially in the THz spectral regime, remain unexplored. Furthermore, exploring the applications of the in-plane HPhPs in optical devices has so far remained elusive. These can be done if one has a vdW crystal with a large enough lateral size while maintaining thicknesses of nanometer scale, so that patterns are made larger than the diffraction limit for these free-space wavelengths and, thus, suitable for far-field spectroscopy.
In this article, we demonstrate far-field excitation and far-field characterization of HPhPs in an In-Plane Hyperbolic Polariton Tuner, which is formed by patterning one-dimensional (1D) ribbon array directly onto the semiconducting HPhP vdW α-MoO 3 flake with a centimeter lateral size while maintaining thicknesses of 100~200 nm. The THz and LWIR photons from far-field illuminating onto the tuner will strongly couple with the phonons and give rise to polaritons with in-plane hyperbolicities. This makes an in-plane tuner an actual device that acts not only with functions of grating but also as a polarizer and notch filter in the LWIR and THz regimes, and have distinctive features including high-Q (300) resonance and extinction ratios up to 6.5 dB at a deep sub-wavelength thickness of 200 nm. Moreover, the polariton resonance frequency, i.e., the operation frequency of the polarizers and notch filters, can be highly tuned by varying the period and the skew angle of the ribbon array.

Fabrication of tuners and far-field HPhPs excitation
An in-plane tuner will have a structured surface written with desirable patterns. In this study, the in-plane tuner consists of simple onedimensional periodic ribbon patterns (1D-PRPs) directly formed on a vdW α-MoO 3 flake using electron-beam lithography (EBL). They have widths (w) and skew angles (θ), which is defined by the angle between the long-axis ribbon direction and [001] crystallographic axis of vdW α-MoO 3 (Fig. 1a). The ribbon period (Λ) is set as 2w. Both of w and Λ are much smaller than the excitation wavelength. As such, we synthesized a 120 nm-thick α-MoO 3 crystal with a largest lateral size larger than 1 cm ( Fig. 1b and Supplementary Fig. 1, and see Methods for details), which guarantees the fabrication, characterization, and comparison of different tuners on the same flake ( Fig. 1c and Supplementary Fig. 2). The basal plane of the as-grown α-MoO 3 is (010) plane, with the two orthogonal directions corresponding to [100] and [001] crystallographic axes, respectively [8][9][10] . In our study, these two axes are defined as the x-and y-axes, respectively (Fig. 1a), which are identified experimentally using micro-Raman spectroscopy ( Supplementary  Fig. 1b-d). A broadband THz and LWIR light illuminates the tuner and the optical responses at the far-field were measured using a polarized Fourier transform infrared (FTIR) micro-spectroscopy (Fig. 1a) (see Methods for details). It should be noted that usually, three main techniques are employed for determining the broadband polaritonic properties of 2D crystals, including the FTIR 36,47,48 , electron energy loss spectroscopy (EELS) 49 , and infrared nanoscopy 4,9,28 . In comparison with the latter two techniques, which usually require complex instrumentation, harsh sample preparation, and time consumption, the farfield polariton characteristics of the 1D-PRPs can be readily measured in a common broadband FTIR spectrometer at ambient conditions, with low time-consuming, a high collection efficiency, and over a large sample area.
In vdW α-MoO 3 crystal, in the spectral regimes 230−400 cm −1 (THz) and 545−1010 cm −1 (LWIR), there are a series of Reststrahlen bands where the Re(ε) along one of the three crystallographic axes, i.e., [100], [001], and [010], is negative ( Supplementary Fig. 3), while at least one is positive. This makes α-MoO 3 a natural hyperbolic medium capable of supporting HPhPs. We first characterized the far-field reflection of the homogeneous pristine α-MoO 3 , which is calculated as R/R 0 − 1, with R and R 0 the reflectance of the light from the surfaces of the sample and bare substrate (see Methods). For incident light polarized along the [001] and [100] directions, the reflectance spectra show distinct peaks at 550 and 820 cm −1 (Fig. 1d), which correspond to IR-active TO phonon modes along [001] ðω 001 TO Þ and [100] axes ðω 100 TO Þ, respectively. Both of the spectra exhibit small valleys at 1004 cm −1 . These valleys are very close to the frequency of LO phonon mode along the [010] axis (z-axis, ω 010 LO ) where the permittivity diminishes. For a 120 nm-thick α-MoO 3 flake, leaky modes (Berreman modes) can be excited near this epsilon-near-zero (ENZ) region and then give rise to the two valleys on the reflectance spectra 50,51 . These narrow leaky modes can further interfere with the broad reflection background and generate asymmetric Fano lineshapes (see Supplementary Note 1 and Supplementary Fig. 4). However, no evident spectral peaks corresponding to HPhPs are observed in the rest of the spectral range. This is due to the large wavevector mismatch between free-space photons (k 0 ) and polaritons (q PhPs ), which prevents the coupling of electromagnetic fields to HPhPs.
A tuner comprised of 1D-PRPs ( Fig. 1c and Supplementary Fig. 2) is able to overcome the large momentum mismatch and excite the HPhPs 52,53 . The incident waves will be scattered by the sharp ribbon edges into evanescent waves with large momenta, whereby HPhPs propagating transverse to the ribbons are excited. Fabry−Pérot resonances (FPRs) can then be formed upon the multiple polariton reflections from the ribbon edges. Simultaneously, the 1D-PRPs can also diffract the incident light into guided waves propagating perpendicular to the ribbon long axis, whose wavevectors are much larger than the free-space waves 54 . These guided waves can then couple with and transfer energy to the polariton FPRs (see Supplementary Note 2 and Supplementary Figs. 5, 6 for more discussion on the excitation of HPhPs by the 1D-PRP). The polariton energy will be dissipated by the lattice vibrations or radiated back to the free space, as manifested by resonance peaks and valleys in the corresponding reflectance and transmission spectra, respectively.
For a typical 1D-PRP with w = 800 nm and orientated along the [001] axis (sample with skew angle θ = 0°shown in Fig. 1c), HPhPs in Reststrahlen Band 2 (820 to 972 cm −1 , where Re(ε x ) <0 and Re(ε y ), Re(ε z ) >0) will be excited upon illumination polarized along the [100] axis. This will lead to a strong reflectance peak at 874 cm −1 , as shown in Fig. 1e. The small bump at 820 cm −1 originates from the intrinsic ω 100 TO . No resonance peaks are observed when the polarization is switched to [001] direction. In contrast, for 1D-PRP orientated along the [100] axis (sample with θ = 90°shown in Fig. 1c), a clear peak at 689 cm −1 is identified (Fig. 1f), suggesting the launching of HPhPs in Reststrahlen Band 1 (545 to 851 cm −1 , where Re(ε y ) <0 and Re(ε x ), Re(ε z ) >0). For polarization along the [100] axis, only the ω 100 TO peak appears. The excitation of HPhPs in both Bands 1 and 2 using these two types of 1D-PRPs can be further confirmed by simulating the near-field  distributions they support. The simulated field distributions at ω = 874 and 689 cm −1 clearly reveal that polaritonic rays with zig−zag shape patterns propagate inside the ribbons ( Supplementary Fig. 7a, b and Supplementary Note 3). These are typical fingerprints of HPhP waveguide modes. These modes are bulk modes with electromagnetic fields confined inside the body of the flakes, such as those observed in hBN nanostructures 55 and biaxial α-MoO 3 flakes 56 .
Notably, in the two 1D-PRPs two valleys with narrow linewidths appear around 1000 cm −1 for both polarization conditions. These valleys are spectrally close, but occur with different amplitudes, with the spectral valleys deeper when the incident light is polarized perpendicular to the ribbon. In addition, as discussed below, the deeper valleys shift when w and θ change. Therefore, they are ascribed to HPhP resonances in Reststrahlen Band 3 (958 to 1010 cm −1 , where Re(ε z ) <0 and Re(ε x ), Re(ε y ) >0). This is also corroborated by the corresponding near-field distributions showing polaritonic rays with zig−zag shapes ( Supplementary Fig. 7c, d). Due to their narrow linewidths, Fano interference will occur between the background reflectance and the HPhP resonances, giving rise to these valley features (see Supplementary Note 1 and Supplementary Fig. 4). The shallower valleys are contributed by the aforementioned ENZ condition that occurs near the ω 010 LO (Fig. 1d). These results are consistent with polariton propagation in the basal plane of an α-MoO 3 flake: the HPhPs in Band 1 and 2 are of in-plane hyperbolicity, which cannot propagate along [100] and [001] directions, respectively 10 . In contrast, the dispersion of HPhPs in Band 3 are elliptical, which therefore allows them to propagate along both of the two orthogonal crystallographic directions 10 .

Tuning the HPhPs with tuners of different ribbon widths
The above results clearly prove that the HPhPs supported by α-MoO 3 , which previously were only accessed using near-field nanoimaging [8][9][10]30,31 , are able to be excited from the far-field using the 1D-PRPs. To demonstrate the tunability of HPhPs in this material, we measured the reflectance spectra of 1D-PRPs with different w ranging from 100 to 2000 nm. All spectra were collected with the incident light polarized perpendicular to the ribbon long axis. For 1D-PRPs parallel to the [100] direction, both of the HPhP resonances corresponding to Band 1 (peaks) and 3 (valleys) can be excited (Fig. 2a). Notably, when w increases from 600 to 2000 nm, the resonance in Band 1 clearly redshifts from 746 to 613 cm −1 , whereas a reverse trend appears for the modes in Band 3, where the resonances blueshift from 998 to 1005 cm −1 . Similar spectral evolution with changing w can be observed for 1D-PRPs along the [001] direction (Fig. 2b), where PhP resonances in Band 2 and 3 are excited. Additionally, the variation of resonance frequency in Band 3 with w is different for these two ribbon orientations ( Supplementary Fig. 8).
The evolution of the HPhP resonances with changing w can be understood by considering that the excitation of HPhPs are originated from the synergy between guided waves of the array and FPRs in an individual ribbon: the scattering of light at the ribbon edges excite the polariton FPR, while the guided waves further couple with and transfer energy to the polariton waves (see Supplementary Note 2 and Supplementary Fig. 5, 6). The conditions for the occurrence of the FPRs and guided waves are q PhPs w = ± mπ and q PhPs Λ = ±n2π, respectively, with m, n = 1, 2, 3, 45,57 . Because in our study, the Λ is deliberately set as 2w, these two equations are the same. Therefore, each w corresponds to an in-plane polariton wavevector of q PhPs = ∓ mπ w . When w changes the resonance peak will scale according to the in-plane polariton dispersion relations, ω(q PhPs ). This can be readily seen by calculating the ω(q PhPs ) along the [100] and [001] directions, which is visualized as the 2D false color plot of the imaginary part of the complex reflectivity Imr p (q PhPs , ω) (see Supplementary Note 4 for details on the calculation of Imr p ). The PhP resonance frequencies obtained from the spectra shown in Fig. 2a, b are then overlaid onto the same plot by using q PhPs = mπ/w. Excellent agreement is obtained between the experimental measurements and calculated lowest-order (m = 1) HPhP branches for all three bands (Fig. 2c). Such an agreement further validates that the highly anisotropic HPhPs in the α-MoO 3 flake are directly excited from the far-field with the help of the in-plane hyperbolic polariton tuners. It is noted that the reflection of polariton waves by the ribbon edges may induce phase shifts, which can violate the condition for the FPRs by a phase of 2Φ 58 . A previous theoretical result showed that in monolayer graphene, the Φ for plasmon polaritons is0 .75π 58 . In comparison with the plasmon polariton in graphene, the propagation and reflection of the in-plane HPhPs are rather complicated, making the phase shift difficult to be predicted. In our analyses, the Φ is taken as π according to a very recent study 59 . The good agreement between the experimental measurements and calculated results further validates our setting.
The far-field reflectance spectra allow for evaluating the Q-factor of the HPhP resonance, which is defined as Q = ω 0 Γ , with ω 0 and Γ, the frequency and linewidth of a specific resonance 7 (see details in Supplementary Note 1 for extraction of the Q-factors). For HPhP resonances in Band 1 and 3, their Q-factors respectively decrease and increase monotonically against ω 0 (Fig. 2d). This is because when ω 0 increases, the PhP resonance in Band 1 shifts closer to the ω 100 TO (820 cm −1 ), while the resonance in Band 3 shifts away from the TO phonon mode along [010] axis at 958 cm −1 ðω 010 TO Þ 10,36 . Thus the polariton dissipation by lattice absorption will be strengthened (suppressed) for Band 1 (Band 3), giving rise to a larger (smaller) Γ. The non-monotonic behavior of the Q-factor for the HPhP resonance observed in Band 2 can be understood by considering that, in addition to ω 100 TO , there is another LO phonon mode along the [100] axis at 972 cm −1 ðω 100 LO Þ 10,36 . Leaky-mode absorption induced by ENZ condition also occurs near the ω 100 LO . Therefore, the Q-factor first increases as ω 0 is farther from the ω 100 TO , and then decreases gradually approaching the ω 100 LO . It is noted that the Q-factors observed in Band 1, 2, and 3 are 15-25, 25-100, and 200-300, respectively. Most of these values are higher than those observed in graphene nano-gratings with similar resonance frequencies 41,42,47 . In particular, the Q-factor of the Band 3 resonance can be as high as 300, which is on par with the highest observed in hBN nanoresonators (360) via the same far-field technique 60 . Such high Qfactors, coupled with the small modal volumes and footprint of the 1D-PRPs, indicate that the α-MoO 3 tuners offer important application potential for high-efficiency compact photonic devices and components, as demonstrated below.
The far-field spectra also allow for extracting the polariton lifetimes, which span from 0.2 to 3.0 ps in the three Reststrahlen bands (Supplementary Note 5, Supplementary Fig. 9, and Supplementary Table 1). It is noted that the fabrication of 1D-PRPs can lead to damage to the ribbon edges, which can reduce the polariton lifetimes and Q factors. Usually, this issue is inevitable during the patterning of the vdW crystal for various characterizations and device applications. To evaluate the additional polariton damping induced by the patterning processes, we performed near-field measurements on the same α-MoO 3 flake in the unpatterned region and extracted the intrinsic polariton lifetimes (Supplementary Note 5, Supplementary Fig. 9, and Supplementary Table 1). In comparison with the unpatterned α-MoO 3 , the phonon polariton lifetime in the ribbon arrays are reduced by 20−57%. The reduction in the lifetime of polariton modes in α-MoO 3 ribbons is attributed to defects or impurities at the rough edges of the ribbons introduced during the fabrication process. These imperfections can create additional scattering centers for HPhPs. Furthermore, they will also lead to the broadening of the pristine phonon modes. Raman spectroscopic characterizations show that, compared to the unpatterned region, the overall linewidths of the phonon modes in the ribbons increased by 4.6-14.1% (Supplementary Note 6, Supplementary Fig. 10, and Supplementary Table 2). This broadening shortens the lifetime of the associated polariton modes as well. Despite the lifetime reduction, the 1D-PRPs still exhibit high Q factors upto 300. Although such a value is smaller than that of resonances sustained by an unpatterned and naturally grown α-MoO 3 ribbon 61 , it can be further improved by optimizing the processing parameters.

Tuning the HPhPs with ribbon arrays of different skew angles
The in-plane HPhP dispersions of the α-MoO 3 are highly anisotropic 9,10 . This offers unique tunability of the HPhP resonances by changing the 1D-PRP orientations, which cannot be realized in vdw crystals with isotropic in-plane dispersions such as hBN and graphene. As such, 1D-PRPs with fixed w (480 nm), but different θ were fabricated (Figs. 3a,  1c). Each pattern can provide polariton momenta of q PhPs = π/w, with the direction perpendicular to the ribbon long axis. In this way, when the ribbon is rotated away from the [001] axis, HPhPs with wavevectors of different orientations within the basal plane can be excited, giving rise to HPhP resonances that are strongly dependent on the θ. Specifically, HPhP resonances in Band 2 and 3 can be observed for θ = 0°, and all of the HPhP resonances in Band 1, 2, and 3 appear when θ is increased (Fig. 3b). For θ = 90°, the resonance in Band 2 merges with the ω 100 TO . In addition, the resonances in Band 1 and 2 are highly sensitive to θ, while that in Band 3 shifts slowly against θ (Fig. 3b). This can be ascribed to the distinct in-plane polariton dispersions of the three bands. Specifically, the in-plane isofrequency contour (IFC) of HPhPs in Band 3 is an ellipse, where the polariton dispersion differs moderately along different directions in the basal plane. However, the IFCs in Band 1 and 2 are hyperbola, making their polariton dispersions highly dispersive with θ. These angle-dependent behaviors can be seen more clearly by plotting the dependence of resonance frequencies in the three bands on the skew angle (Fig. 3c), which agrees well with the calculated Imr p (q PhPs , ω).  The HPhP dispersion relations at each skew angle can be obtained by measuring the reflectance spectra from 1D-PRPs with different w (q PhPs ) at a specific θ ( Supplementary Fig. 12), whereby the in-plane polariton IFCs at different energies can be re-constructed and visualized. Clearly, the HPhP resonances in Band 1 and 2 depict IFCs of open hyperbolic shapes (Fig. 3d, e), while those in Band 3 correspond to an IFC of a closed ellipse (Fig. 3f). Moreover, at higher frequencies in Band 1 (Band 2), the opening-angles of the hyperbolic sectors become smaller and the hyperbola bends toward the [001] ([100]) direction (Fig. 3d, e). All the experimental points can be fit well by the calculated Imr p (pseudo-colored plots shown in Fig. 3d-f and Supplementary Fig. 12). These results provide further direct evidence for far-field excitation and modulation of the hyperbolic HPhPs in the α-MoO 3 flake.
High-Q resonances can also be induced in 1D-PRPs made out of an in-plane isotropic polaritonic vdW crystal 44 . The frequency of such resonances can also be tuned by changing the ribbon width, but the α-MoO 3 1D-PRPs proposed in the current study is unique. The in-plane hyperbolicity of α-MoO 3 makes it possible to tune the resonance frequency of the 1D-PRPs by rotating the ribbons while fixing their widths. This feature can greatly simplify the fabrication processes of arrays with different resonance frequencies, and small additional damping will be introduced because the ribbon width is unchanged. Moreover, for a polariton mode with a wavevector approaching the asymptote of the hyperbolic IFC, the polariton momentum will become remarkably high. This will generate much stronger electromagnetic field confinement than those with wavevectors away from the asymptote. These   modes can significantly enhance light−matter interactions at the nanoscale and lead to various applications, such as enhanced light emission, ultrasensitive biosensing, and nonlinear optical signal generations. For these applications, usually, a fixed operation frequency is preferred. In the 1D-PRPs made out of α-MoO 3 , tuning the orientation and width of the ribbons can both tailor the resonance frequency. Therefore, it is possible to induce a high-momentum polariton mode while fixing its resonance frequency by simultaneously orientating the ribbon long axis along the asymptote of the hyperbolic IFC and tuning the ribbon width. This feature can open an avenue for the applications just mentioned.

Far-field excitation of tunable THz HPhPs with ribbon arrays
The α-MoO 3 flake can also support nanoscale-confined HPhPs in the THz domain from 260-400 cm -1 (8-12 THz), which, however, have only been probed using the near-field nano-imaging technique 30 . To demonstrate the far-field excitation and tuning of the THz HPhPs, we fabricated 1D-PRPs of different w and θ and characterized their spectral responses. Due to the relatively low signal-to-noise ratio of the bolometer in the THz domain, transmission spectra were recorded, which is defined as T/T 0 , with T and T 0 the transmittance of the light through the sample and bare substrate (Methods  (Fig. 4a). The tunability of these two resonances is clearly demonstrated by their redshifting behaviors with increasing w. The polariton dispersion relations were then obtained by extracting the resonance frequencies at different w. The analytical dispersions were derived using the dielectric tensor reported in ref. 30 and are in good agreement with the experiment measurements ( Fig. 4b and Supplementary Note 7). The Q-factors of the two types of resonances can be evaluated according to the transmittance spectra shown in Fig. 4a, which both increase with increasing the resonance frequency (Fig. 4d). The available Q-factors are in the range of 15-25 (HPhP [001] ) and  ) respectively, which surpass that of graphene plasmon resonance in the THz regime 40 . The HPhP resonances are also strongly dependent on the skew angle of the 1D tuner patterns. By sweeping the θ from 0°to 90°( ribbon long axis rotating from [001] to [100] direction), the HPhP [001] and HPhP [100] resonances shift monotonically to higher and lower frequencies, respectively (Fig. 4c). Moreover, with the experimental data shown in Fig. 4b, c, the in-plane IFC of the HPhP [001] can readily be drawn to exhibit a clear hyperbola opening towards the [001] axis ( Supplementary Fig. 13). These results corroborate the excitation of HPhP resonances in THz Reststrahlen Band 1 (HPhP [001] ) and 3 (HPhP [100] ), which are consistent with the previous nano-imaging results 30 .

Tunable LWIR and THz polarization notch filters
A PNF is a unique optical component that combines a polarizer and a narrow band-rejection filter together into a single component, which is able to block a monochromic laser with a given linear polarization, while passing light of all polarization states at wavelengths adjacent to the laser line. Such a filter has broad application prospects in laser spectroscopy and optical communications, but the commercial products are rare, and especially, there is no commercial PNF in the LWIR and THz ranges. The salient high-Q (Figs. 2d, 4d) and polarizationsensitive (Fig. 1e, f) α-MoO 3 tuners established above provide opportunities for developing tunable PNFs 62 . As such, we constructed PNFs using 1D-PRPs with long axes along [001] and [100] axes, respectively. For typical 1D-PRPs with w = 1000 nm, their extinction spectra in LWIR (Fig. 5a) and THz (Fig. 5b) regimes are strongly polarizationdependent. Specifically, the HPhP resonances only appear when the excitation polarization is perpendicular to the ribbon's long axis. This can be seen more clearly by plotting the extinction at the corresponding resonance peaks, i.e., 650 cm −1 /270 cm −1 and 869 cm −1 / 362 cm −1 , against light polarization for ribbons parallel to [100] and [001] directions, respectively (Fig. 5c, d). The performance of a PNF for blocking a monochromatic light source can be evaluated by two parameters: the polarization extinction ratio and bandwidth. Specifically, the polarization extinction ratio is defined as 10 log(T 0 /T), with T 0 and T the transmittance of light polarized along and perpendicular to the ribbon's long axis at the resonance frequency. The bandwidth can be calculated according to the full width at half maximum (FWHM) of a specific resonance. Accordingly, for the PNFs with resonances at 650, 869, 270, and 362 cm −1 , corresponding extinction ratios (bandwidth) are 3.4 dB (41.5 cm −1 ), 5.5 dB (26.5 cm −1 ), 4.5 dB (23.4 cm −1 ), and 3.7 dB (24.5 cm −1 ), respectively. In particular, the bandwidths of the PNFs are comparable to many of the commercial narrow-band-pass filters in similar spectral regimes, which are usually polarization insensitive, (Supplementary Fig. 14), while the thicknesses of the PNFs (200 nm) are much smaller than those of the commercial components (~1 mm). More importantly, using 1D-PRPs with different w, the operational frequency, polarization extinction ratio, and FWHM of the PNF can be engineered continuously ( Supplementary Figs.15, 16). The maximum peak extinction can reach 6.5 dB and the smallest FWHM can be as narrow as 17 cm −1 .

Discussion
We have successfully demonstrated direct far-field excitation and characterization of the tunable LWIR and THz HPhPs in biaxial vdW α-MoO 3 patterned into simple 1D-PRPs. The 1D-PRPs can act as polariton tuners that are sensitive to the excitation polarization and with light extinction ratios up to 6.5 dB and high Q-factors up to 300. Such a compositional set of output functions are tunable and strongly dependent upon the in-plane hyperbolic phonon polaritons in the α-MoO 3 . It is noted that in comparison with the recently reported patterned vdW WTe 2 flakes with tunable in-plane hyperbolic plasmon polaritons at cryogenic conditions 45 , the polariton tuners proposed in the current study may be more favorable for practical applications because of their room-temperature operation, broader spectral range, and much higher quality factors of the resonance modes.
On the prospects, these in-plane hyperbolic polariton tuners can be used in optical circuitry, instruments and even modern information systems. The in-plane hyperbolic polariton tuner also opens up avenues for a variety of practical photonic and optoelectronic applications besides the PNFs shown in Fig. 5. For example, by engineering the tuner to spectrally overlap the HPhP resonance with a specific vibration or rotation transition of a molecule, strong interactions between the tuner and molecule can be induced 44,63 , which can significantly enhance the molecular absorption or emission and give rise to various ultrasensitive bio-sensing techniques. The tuners with high-Q and tunable HPhP resonances can also be employed to regulate the blackbody emission 62 , whereby narrow-band, polarized, and tunable thermal emission can be achieved [64][65][66] . Moreover, tunable and highperformance LWIR and THz photodetectors can also be envisioned by taking advantage of the strong light field localizations ( Supplementary  Fig. 7) and semiconductor nature of the α-MoO 3 19 ; this type of devices may give unique functions necessary for future communication and radar applications.
For fundamental research, the far-field excitation methodology can complement near-field nano-imaging techniques and make the characterizations of PhPs of materials more precisely, especially for those with in-plane hyperbolicity. For example, in nano-imaging measurement, because the polariton waves are launched by an antenna (the scanning tip or antenna fabricated onto the sample surface), in principle, these HPhPs can propagate along different directions determined by the hyperbolic IFC [8][9][10] . Therefore, the measurements of the polariton wavelength and dispersion relation from the polariton interference fringes can be disturbed by these various polariton waves. On the other hand, with far-field excitation, the polariton wavevector is determined by the tuner's structural parameter and orientation. For a specific pattern, only one HPhP mode can be excited, which allows for a more accurate characterization of the intrinsic HPhP properties. Additionally, it is possible to characterize the HPhPs using a broadband light source covering a broad THz spectral range. This can unveil more complete polaritonic properties for broadening and deepening our understanding of the THz HPhPs, in particular for two-dimensional  atomic crystals (e.g., the in-plane IFCs in the THz domain as shown in Supplementary Fig. 13), which is now limited by the discrete laser lines used in most nano-imaging measurements 30 . It is noted that in the current study, as a demonstration of principle, we only employ the simplest form of patterns, and only demonstrate with one type of material. In fact, patterns can be chosen depending on desirable applications, and also materials as long as supporting HPhPs in the THz and LWIR regimes 6,67 . The tuner allows us to modulate the wavefront of the incident light and control their power flow in an engineered space [21][22][23][24][25][26] , which therefore enables a variety of interesting applications, such as negative refraction, holography, metalens, polarization conversion, and even topological PhPs and exciton-polaritons with robust beam steering functionalities 68,69 , in visible and near-infrared ranges to be expanded into the THz and LWIR spectral regimes.
When preparing the revised manuscript, we became aware of a recent work 59 on a similar topic to our current study.

Methods
Synthesis of large-area vdW α-MoO 3 flake vdW α-MoO 3 thin flakes with a centimeter lateral size were prepared using a modified thermal physical vapor deposition method. Specifically, an alumina crucible filled with 0.1-g α-MoO 3 powders was placed at the center of a quartz tube as the evaporation source. Another crucible covered with a silicon substrate of 1 cm × 1 cm was placed 20 cm away from the source. The source was then heated up to 780°C  [001] and held at that temperature for 2 h. The α-MoO 3 powders were sublimated and crystallized onto the silicon wafer. Afterward, the quartz tube was cooled down to room temperature naturally. The large-area α-MoO 3 flakes can be found on the silicon wafer.

Fabrication of the vdW α-MoO 3 1D periodic tuner patterns
The as-grown α-MoO 3 flakes were transferred to a pristine (highly resistive) silicon substrate covered with a 300-nm oxide layer. Afterward, a selected flake was patterned into 1D-PRPs with different w and θ using a combination of electron beam lithography (EBL: EBPG5000+, Netherlands) and reactive ion etching (RIE: 50 W for 10 min). For the EBL processing, a 400-nm layer of Poly(methylmethacrylate) (PMMA) photoresist was used. For the RIE etching, a mixture of O 2 (12 vol.%), Ar (30 vol.%), and CHF 3 (58 vol.%) was employed. The etching was conducted at 50 W for 10 min. To guarantee good signal-to-noise ratios of the spectral characterizations, the areas of the patterns were set as 50 μm × 50 μm and 300 μm × 300 μm for LWIR and THz regimes, respectively.
Reflectance and transmittance spectral characterizations LWIR and THz spectral characterizations were performed using a Bruker FTIR spectrometer (Vertex 70 v) integrated with a Hyperion 3000 microscope and a mercury cadmium telluride (HgCdTe) photoconductor (for the measurement of LWIR spectra) or a liquid-helium-cooled silicon bolometer (or the measurements of THz spectra) as the detector. A broadband black-body light source covering the LWIR and THz spectral regimes was employed as the incident light. A 15× reflective Schwarzschild objective was utilized to focus the incident light onto the tuner patterns with a spot size of~500 μm. The reflected or transmitted light were collected from an area of 50 μm × 50 μm and 300 μm × 300 μm for LWIR and THz regimes, respectively, with the help of an iris diaphragm. For the polarization-dependent measurements, a polarizer was used to control the polarization of the incident light. For the reflectance and transmittance spectra of the tuner, a bare silicon substrate is used as reference for normalization.
To characterize the PNFs in LWIR and THz regimes, a linear polarizer was placed before the detector to determine the polarization state of light transmitting through the tuners upon an unpolarized illumination.

Data availability
Relevant data supporting the key findings of this study are available within the article and the Supplementary Information file. All raw data generated during the current study are provided in the Source Data file or available from the corresponding authors upon request. Source data are provided with this paper.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/.